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Improving the precision of light quark mass determinations 

Christian Sturm 
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We discuss the concepts and the framework of the renormalization procedure in regularization- 
invariant momentum subtraction schemes. These schemes are used in the context of lattice sim- 
ulations for the determination of physical quantities like light quark masses. We focus on the 
renormalization procedure with a symmetric subtraction point of a quark mass and discuss its con- 
version from the RI/SMOM scheme to the MS scheme. A symmetric subtraction point allows for a 
lattice calculation with a reduced contamination from infrared effects. The perturbative series for 
the resulting matching factor at next-to-leading order is better behaved for the symmetric momen- 
tum configuration than for the exceptional one which can decrease the uncertainty in light quark 
mass determinations. 



I. INTRODUCTION 

In continuum perturbation theory dimensional reg- 
ularization [l[ is for many phenomenological applica- 
tions the most common regulator for divergent loop- 
integrals. Physical observables computed to higher 
orders are often conveniently rcnormalizcd in the min- 
imal subtraction (MS) scheme [2|, |3j; thus the renor- 
malization constants are well known to high order in 
perturbation theory. 

However, the MS scheme is fixed to dimensional reg- 
ularization and is usually not directly applicable in 
lattice calculations. In the context of lattice simula- 
tions in combination with non-perturbative renormal- 
ization, regularization-invariant momentum subtrac- 
tion(MOM) schemes are used 4]. In regularization- 
invariant (RI) schemes the renormalization constants 
are often fixed by demanding that the renormalized 
quantity is equal to its lowest order value for a special 
momentum configuration(subtraction point). In order 
to convert results computed in a RI/MOM scheme to 
the MS scheme matching factors are needed which 
translate an observable renormalized in the given 
RI/MOM scheme into the MS scheme. These match- 
ing factors can be computed perturbatively. For ex- 
ample, the computation of the mass conversion fac- 
tor C^^ MOM , which turns a quark mass renormal- 
ized in the RI/MOM scheme into the MS scheme, 
or the conversion factor (j^- l / MOM ; wn ich performs 
the corresponding conversion of the quark fields, are 
both known up to three-loop order inperturbative 
Quantum Chromodynamics(QCD) [1, I n the 

Rl'/MOM scheme which is also used in lattice simu- 
lations these matching factors are also known to the 
same order 0, 0| . 

Apart from quark masses, also the strong coupling 
constant a s has been considered in momentum sub- 
traction schemes @, ©, [H El, El. Taking into 
account vertex diagrams, one has a larger choice of 
defining the subtraction point at which the renor- 
malization constants are fixed through different mo- 
mentum configurations. Two examples are the excep- 
tional^!' asymmetric) subtraction point with the ex- 



ternal minkowskian momenta p\ = pi, = —[J. 2 , P3 = 
or the non-exceptional(or symmetric) momentum con- 
figuration with p\ = p\ = p\ = ~M 2 - The symbol fi 
denotes here the renormalization scale. 



In Ref . [14| the renormalization of quark bilinear op- 
erators for a symmetric momentum subtraction point 
has been studied in the RI/SMOM scheme for the 
vector, axial- vector, scalar, pseudo-scalar and tensor 
operators, where the S in SMOM stands for sym- 
metric. The corresponding renormalization constants 
Z are related among each other through the Ward- 
Takahashi identities. These relations have been used 
to perform mass renormalization with a symmetric 
subtraction point. 

Light up-, down- and strange-quark masses were de- 
termined in Ref. [l5| through lattice simulations in 
the RI/MOM scheme and subsequently converted to 
the MS scheme. However, the perturbative expan- 
sion of the matching factor C^^ MOM in QCD shows 
poor convergence and leads currently to a significant 
contribution to the error on these quark masses. The 
uncertainty from the renormalization procedure of the 
quark masses in Ref. [15[ amounts to around 60% of 
the total error. In order to renormalize the quark 
masses in the lattice simulation, the renormalization 
constants need to be determined in a lattice calcula- 
tion [l6| . As shown in Ref. , a symmetric subtrac- 
tion point involves non-exceptional momenta which 
imply a lattice simulation with suppressed contami- 
nation from infrared effects, whereas the asymmetric 
subtraction point is less suited, since effects of chiral 
symmetry breaking vanish only slowly for large exter- 
nal momenta as a result of Weinberg's theorem [l7j |. 

This proceedings contribution is organized as fol- 
lows: in Section|TT]we introduced the used conventions 
and notations and discuss generalities. In Section [TTT1 
we summarize the results for the mass renormaliza- 
tion with a symmetric subtraction point, addressing 
the matching factors and the two-loop anomalous di- 
mensions. Finally we close with the conclusions in 
Section [TV] 
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II. GENERALITIES AND NOTATION 

Bare and renormalized quantities are related 
through the renormalization constants Z, like the 
renormalization constant Z m of the quark mass m and 
the renormalization constant Z q of the fermion field 5* 
with the properties 



tub. = Z m m,B and 



(1) 



The subscripts R and b denote renormalized and bare 
objects. These renormalization constants can be fixed 
through renormalization conditions imposed on self- 
energy type diagrams as shown in Fig. Q] given by 



iS(p) = dxe tpx {T[^(x)^{Q)]) 



' — m + ie — 



(2) 
(3) 



where incorporates the higher order QCD cor- 

rections. 




(a) 



(b) 



FIG. 1: QCD corrections up to one-loop order to the quark 
propagator. 

In perturbation theory the conditions for the self- 
energies in the last line of Table U fix the renormal- 
ization constants in the RI/MOM scheme. 

The vector and axial-vector Ward-Takahashi iden- 
tities relate the self-energies and amputated Green's 
functions 



A 6 = S~ 1 ( P2 )G 6 S- 1 (p 1 ) 



(4) 



of quark bilinear operators O = uT^d with the in- 
sertion of, in general, the vector(ry = 7^), axial- 
vector(r j 4 = 7^75), scalar(rs = 1) or pseudo- 
scalar(rp = 175) operators. Diagrams of these 
Green's functions are shown in Fig. [5] 

For each of the operators a renormalization constant 
Zq needs to be introduced with Op = Z Q Op- In the 
MS, RI/MOM and RI'/MOM schemes the renormal- 
ization constants are related among each other with 
Zy = 1 = Za, Zp = 1 jZ^ and Z§ = Zp; more details 
can be found in e.g. Refs. [TalHI. Some of these prop- 
erties hold non-perturbatively while others are proven 
only in continuum perturbation theory. Through the 
Ward-Takahashi identities the renormalization con- 
stants in the RI/MOM scheme can also be determined 




3- Hi H 2 



FIG. 2: Diagrams which contribute to the computation of 
the Green's function Gq at leading order (first line) and 
next-to-leading order(second line) in perturbative QCD. 
The circle with the cross denotes the inserted operator O. 
Solid lines indicate fermions and spiral lines gluons. 



by the conditions on the amputated Green's functions 
as given in the last row of TableQ] In this case the sub- 



traction point is exceptional, pf 



Pi 



-AT, p 3 = 0, 



without any momentum transfer leaving the operator. 



In Ref. 1J| a symmetric subtraction point has been 



studied and its applicability demonstrated through 
the RI/SMOM and RI/SMOM^ schemes. The con- 
ditions for these new schemes are also summarized in 
Table U and the proofs for the above relations among 
the renormalization constants has been performed in 
Ref. [3. A non-perturbative test of the RI/SMOM 
scheme for the quark mass renormalization can be 
found in Ref. [H- 

Also the tensor operator, with Ft = cr^ = 
|[7' i ,7 l l_has been studied in Refs. 0, M, H| and 
in Ref. :14j] for a symmetric subtraction point. 



III. RESULTS 



A. Matching factors 



The matching factors which convert a quark mass 
from the RI/SMOM or RI/SMOM^ scheme to the 

MS scheme 



7m 



„ ? MS _ ^/SMOM^, RI/SMOM (7m) 

,n R — ,n R \ a ) 

have been computed to one-loop order in perturbative 
QCD in Ref. [14]. The numerically evaluated result 
allows for a direct comparison of the size of the per- 
turbative coefficients in the different schemes. In the 
RI/SMOM scheme the result reads 

c m/smom = i_ ^[0.645518856... 



47T 



-£0.229271492...] + 0(a 2 s ).(6) 
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RI/SMOM 



RI/SMOM 



RI/MOM 



^Tr[g M A^( Pl ,p 2 )^| =1, lim ^Tr [5^^(^2)75 

sym mR^v 

lim ^Tr[A 5 ,fl(pi,P2)D]| = 1, lim ^Tr [Ap,fl( Pl ,p 2 ) 75 ] | = 1, 

1 m rrr-1 I I I 



lim 



777 R 
lil 

rn.fi— *0 r 

lim „T2- 

777fj-»0 



= 1, 



sym 



lim ^ 

77l H ->0 

lim ^Tr[A s ,_ 

777^ — »U 



m R 

W^^R 2 2 = _1 ' 

^ {it [S- 1 ^)] | p2= _ M2 - i Tr [g M A^( Pl , P2 ) 7; 
TrrA^ R (p 1 ,p 2 ) 7 J| = 



i Tr [ A v,h(pi.P2)7m] 
>Tr[A £ 

lim mi? _ 35 | r 



m R ^0 



= 1, 

sym 

«bi,P2)D]| =1, lim ^TrlAp, 



-I I sym J 



1 > JL J I sym J 

l™ n ^Tr [A^ ifl ( Pl ,p 2 ) 757M ] | sy = 



= 1, 



Tr 



l™777 J? -.o T2 ^{Tr[S- 1 (p)]| 



O 2 — — /J,' 

■fi = — yp- 



R (Pl,P2h 5 ]\ sy =h 

= -i, 



lim n s 

m R ->0 

lim ^Tr [As, 

77l H ->0 



+ Tr[ gM7 ^A^]| sj/mj 

51,^2)75] I r = L 

' ' J I sym J 



1 TV 

2 ±L 



lim ^Tr 

m R ->0 



iTr [A^ B (p 1>P2 ) 7M ] ^ = 1, m lhn o ^Tr [A^( Pl , P2 ) 757/1 ] 1, 

Wfa'P^la^ X ' 77^ ^ Tr[AP ^ (Pl ' P2)75] L S!J 77T L 



p 2 = -u- 2 



777 H - 

= -1, 



lim , 

m^— >0 rt 



= 1. 



TABLE I: This table summarizes the renormalization conditions for the RI/SMOM and RI/SMOM^ schemes as defined 
in Ref. [13 ] for the amputated Green's functions of the vector, axial-vector, scalar andpseudo-scalar operators as well 
as for the self-energies. The last row contains the conditions for the RI/MOM scheme [j. The subscript (a)sym stands 
for the restriction to the (a)symmetric momentum configuration and q — pi — p 2 is the momentum transfer leaving the 
operator. 



The corresponding conversion factor for the fermion 
fields cf /SMOM = Zp/Zf /SMOM in the case of the 
RI/SMOM scheme is given by 

c ri/smom = 1 _ £^1.333333333... + (a 2 s ). (7) 

In the RI/SMOM^ scheme the results for the match- 
ing factors read 

RI/SMOM, Q< r 

Cm M = l - 1.978852189... 

47T 

+£0.666666667...] + O (a 2 s ) (8) 



and 



RI/SMOM Q< , 

C q ^ = 1 + — 1.333333333. . 

4-7T 

-£0.437395174. 



O(o2),(9) 



where £ is the gauge parameter and a s the strong 
coupling constant. 

In order to determine the perturbative mass con- 
version factors C m for a symmetric subtraction point 
in the two mentioned schemes, one can consider the 
computation of the amputated Green's function of 
the (pseudo-) scalar operator or the (axial-)vector op- 
erator with suitable projectors. In the case of the 
pseudo-scalar and axial-vector operator a naive anti- 
commuting definition of 75 is used for the treatment 
of 75 in dimensional regularization [J, [23[ which is a 
self-consistent prescription for the flavor non-singlet 
contributions I24l [2511 . 



Comparing these matching factors with the conver- 
sion factor in the RI/MOM scheme 

^ri/mom = i_ ^[5.333333333... 

47T 

+£2.000000000...] + O (a 2 s )(10) 
and the Rl'/MOM schemed!] 

c ri'/mom = 1 _ ^[5.333333333... 

47T 

+£1.333333333...] + O (a 2 s ),(ll) 

one finds in particular in the Landau gauge (£ = 0), 
which is usually adopted in the lattice simulations, 
a significantly smaller one-loop coefficient which sug- 
gests a better behaved perturbative series. If this be- 
havior persists at higher orders it allows for a sub- 
stantial reduction of the systematic uncertainty in the 
light quark mass determinations. 

The matching factors C m in the x = RI/MOM, 
x = Rl'/MOM, x = RI/SMOM and x = RI/SMOM^ 
schemes are gauge dependent. While performing the 
conversion of the quark mass to the MS scheme this 
gauge dependence is compensated by the gauge de- 
pendence of the renormalization constant determined 
in the lattice simulation. The one-loop coefficient 



4tt 



Cm' X {Q of the matching factor C m = 1 
as a function of £ is shown in Fig. [3] for the ex- 
ceptional and non-exceptional momentum configura- 
tions. For almost all gauges the one-loop coefficient 
of the symmetric subtraction point is smaller than in 
the asymmetric case except for a small window with 
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FIG. 3: This figure shows the one-loop coefficients c„ 
of the matching factors Cf n as a function of the gauge 
parameter £ for two exceptional (dashed line(red), x = 
RI'/MOM; dotted line(blue), x = RI/MOM) and two 
non-exceptional (solid line(red), x — RI/SMOM; dashed- 
dotted line(blue), x = RI/SMOM^ ) momentum configu- 
rations. The horizontal lines indicate the size of the one- 
loop coefficient in the Landau gauge(£ = 0). 



£ £ (-5.41532..., -3) comparing the RI'/MOM and 
RI/SMOM schemes and f 6 (-2.74207..., -2.51586...) 
comparing the RI/MOM and RI/SMOM^ schemes. 



B. Two-loop anomalous dimensions 



The mass anomalous dimension 

djogmpj) _ (o)^£_ ( 

7m _ 7 
TT 



dlog {n 2 ) 
1$ {^y + 0{ai) 



i) (^1 
1 \ir 



(12) 



can be used to run the quark mass to different en- 
ergy scales. In the MS scheme the lowest expansion 
coefficients read 



7m — 4° F ' 



(13) 



7 



(D,ms = ± ( 3^ + 97 CfCa _ ^ CFTFnf l (14) 



16 V 2 



6 



10. 

y 



where Cf = 4/3(C^ = 3) denotes the Casimir op- 
erator in the fundamental(adjoint) representation of 
SU(3). The symbol n/ represents the number of ac- 
tive fermions and Tp = 1/2 is the normalization of 
the trace of the SU(3) generators in the fundamental 
representation. 

The mass anomalous dimensions of the RI/SMOM 
and RI/SMOM 7 schemes up to two- loop order are 
given by [lj] 



7 



(0),RI/SMOM _ (0),MS _ 
7m 



(0), RI/SMOM 



In 



',(15) 



(l),RI/SMOM = (1),MS 
/ m J m 



+ /3 (0) C F 0.121034785..., (16) 



(1), RI/SMOM 



7n 



if = ^(1),MS 
/ m 



+ /? (0) C F 0.371034785..., (17) 
with the QCD /3-function defined by 
da s (ii)/TT 







d\og([i 2 ) \ IT 

and the lowest order coefficient 

^ = \(fc A -lTpn f 



/3<°>(^) +0(al) (18) 



(19) 



Similarly one can define the anomalous dimension 
dlogtffl 



7? 



J dlog(pt 2 )' 



(20) 



As shown in Ref. [lj] it is equal in the RI/SMOM 
and RI'/MOM schemes and known up to order o? s in 
Refs. HQ. In the RI/SMOM scheme it is given in 
the Landau gauge by 



2 / 3 
32' 



31 
192 



CfCa 



--C F T F n f )+0(ai) (21) 



and has been determined in Ref. [141 ]. 



IV. CONCLUSIONS 

We have discussed the renormalization of quark 
masses in a regularization-invariant momentum sub- 
traction scheme with a symmetric subtraction point. 
The scheme is useful in the context of lattice simula- 
tions of the light up-, down- and strange-quark masses 
as an intermediate step before the conversion to the 
MS scheme. The renormalization constants can be de- 
termined by the computation of amputated Green's 
functions with the insertion of the scalar, pseudo- 
scalar, vector and axial-vector operators. This allows 
one to derive in continuum perturbation theory the 
two-loop anomalous dimensions and one-loop match- 
ing factors which convert the light quark masses renor- 
malized in the RI/SMOM scheme into the MS scheme. 
The new RI/SMOM scheme decreases chiral symme- 
try breaking as well as other unwanted infrared effects 
in the lattice calculation and the one-loop coefficient 
of the new matching factor is smaller compared to the 
one of the RI/MOM(Rl'/MOM) scheme. A better be- 
haved perturbative series can improve the accuracy of 
light quark mass determinations, if this behavior is 
confirmed at higher orders. 
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